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Abstract 



In this work we apply the Lie group representation method introduced in the 



real time formalism for finite-temperature quantum-field theory, thermofield 

5-H ' 



dynamics, to derive a spinorial density matrix equation. Symmetry prop- 
erties of such equation are analysed, and as a basic result it is shown that 
one solution is the generalised density matrix operator proposed by Heinz, 
to deal with gauge covariant kinetic equations. In the same context, prelim- 
inary aspects of a Lagrangian formalism to derive kinetic equations, as well 
as quantum density matrix equations in curved space-time, are discussed. 
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I. INTRODUCTION 



This paper presents a group symmetry scheme in order to treat the relatisvistic density 
matrix including the gauge covariance. In a previous paper jl] (to be referred from now 
on as i), we have set forth a representation theory for Lie groups, which is based in the 
thermofield dynamics (TFD), a real-time finite-temperature field theory via an operator 
method P-Qj. In I the TFD Hilbert space has been used as the representation space for 
Lie algebras, and so, the equations of motion in quantum TFD have been derived from the 
Galilei and Poincare symmetries [0,0; a classical counterpart of TFD has been identified 
fU?]; and a connection of TFD with the GNS (Gel'fand, Naimark and Segal) construction 
|| has been analysed in particular in association with elements of Hopf algebras |5|,p|-[l2|| . A 
density-matrix equation for the Klein-Gordon field has been introduced and basic aspects 
of the kinetic theory for bosons has been analysed ||. This has indicated the possibility to 
extend this symmetry group formalism for other systems of the quantum kinetic field theory, 
in particular by considering gauge fields and curved space-time. Such a study can bring the 
relativistic transport equations to a more fundamental level of symmetry group, and such 
aspect find interest in the recent developments in high energy physics and in early universe 
theories. 

One remarkable prediction of the quantum chromodynamics is the possibility of a tran- 
sition from a confined hadronic matter to a deconfined quark-gluon plasma (QGP), in an 
ultra-relativistic collision of heavy ions, at an energy density of ~2 GeV/fm 3 and a tem- 
perature of ~ 200MeV. Experimentally, such conditions will be realized soon, and this fact 
makes it essential to elaborate on some aspects, such as the signature of the QGP in the 
pre-equilibrium stage ||13|| . These developments have provided a strong motivation for a re- 



vival of the kinetic theory in order to treat classical or quantum quark-gluon systems [13,14 



It is our contention that the symmetry properties can provide an important tool for such 
developments. Therefore, this paper explores the fundamental consequences of symmetry 
on the kinetic theory under gauge invariance as in the case of a quark-gluon plasma. 
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In order to deal with gauge kinetic-theory models, Heinz fl5|-p2| has proposed a gener- 
alization of the concept of density matrix to take into account explicitly the gauge fields. 
Considering a field (j)(x), Heinz has written the density matrix in the following form 

*(x, y) = e~ yd{x) <j>(x) <g> 4>{x) e yd{x) ' , (1) 

where d(x) is a covariant derivative. Such a theory has been successful in introducing gauge 
invariant Wigner operators in relativistic phase spaces, such that transport equations have 
been derived, though, the collision terms have not been treated in a full and consistent 
way (indeed, difficulties with the collision terms are a characteristic aspect of the present 
state of the art of the gauge kinetic theory). One way to proceed in order to derive a full- 
collision kinetic approach is, first of all, to understand the underlying structure attached to 
the Heinz density operator. One such answer is presented in this paper, where we derive (|]) 
as a solution of a density matrix equation proposed here via the representation theory of 
/ considering the spinorial field. This result suggests a generalization of the Heinz density 
operator (|T]) to deal with quantum fields in curved space-time. 

An approach to derive transport equations for systems treated by a quantum field in 



curved space-time background was proposed by Calzetta, Habib and Hu PB[ and Calzetta 



and Hu ||24j| . In such a method problems related to the notion of particles in curved manifold, 
averages, distribution functions and distance between two points, all of them necessary to 
build a kinetic theory, have been partially overcome through the use of the Riemann normal 
coordinates, in which, in the neighborhood of the origin of coordinates, some modes are 
close to plane waves. So the metric tensor g^ v is expanded as 

where rj^ is a flat-space metric and Rp TVp is the Riemann tensor. On the other hand , 
information about early universes may be obtained by the time evolution of the density 
matrix of the universe, by considering the interaction of the quantum fields and many- 
particles systems with the classical gravitational field. An important problem in this context 
is of the gauge quantum kinetic theory, which we are interested to consider at this stage. 



We proceed in the following way. In section 2 some structural elements of TFD are 
presented in connection with the representation theory developed in /, but rederiving the 
basic method along with our proposals here. In section 3 the density matrix equation is 
introduced via a representation of the Poincare group and the gauge invariance is analysed 
in this context. In section 4 we derive a collisionless Boltzmann equation for the spinorial 
field in an algebraic way, such that the approach proposed here can be better compared with 
others in the literature. A possible Lagrangian formalism for the kinetic theory is discussed 
and some comments are made about the extension to the curved space-time as suggested by 



Calzetta, Habib and Hu p3|,P4|. Final remarks and conclusions are presented in Section 5. 



II. THERMOFIELD DYNAMICS AND REPRESENTATION THEORY FOR LIE 

GROUPS 

The thermofield dynamics is based on two main ingredients: one , the dual (or tilde) 
conjugation rules, introduces a doubling in the dynamical variable; while the other, the Bogo- 
lioubov transformation, introduces thermal effects through a vacuum condensation process. 
The doubling, which is a very attractive algebraic aspect in TFD, is defined in the following 
way. 

Let C = {Ai, i = 1,2, 3, ...n; a, b, c, ... G C} be an associative algebra on a C-field. Tilde 
conjugation is a mapping r : C, denoted as tAt = A, furnishing the following properties 

I: 

(M) = M> (2) 
{cAi + A J ) = c*A i + Aj, (3) 

(4) = (4) 
{Ai) = A, (5) 
[Ai,Aj}=0, (6) 

where [ , ] is the commutator. Usually, in TFD, there are two more rules: \0(0)) = \0(8)) 



and (0(9)\ = (0(6)\, which show that vacua in TFD are tilde-conjugation invariant M. 

Automorphisms in £ can then be introduced through unitary operators, say £/(£), such 
that we require r£7(£)r = t/(£)- With this result, writing £7(£) as = exp(i^A), where 
A is a transformation generator, we have tAt = —A. Therefore, A can be considered as an 
odd polynomial function of A — A, i.e., 

oo 

A = f(A - A) = £ c n (A - Af n+ \ (7) 

n=0 

where the coefficients c n are c-numbers. 

We consider here the simplest case with n = 0, that is, 

A = A - A. (8) 

Therefore, this, Eq.(^), equips the TFD algebra, C, with an hat-isomorphism, satisfying the 
properties: 

(cA + BY= cA + B, c G R, (9) 

{icAY= ic(2A - A), (10) 

(AB)~= AB + AB - AB, (11) 

[A,B]-=[A,B] + [A,B]-[A,B], (12) 

[A, B) + ~= [A, B} + + [A,B) + - [A,B) + , (13) 

A = 2A. (14) 

The proof follows from the definition of tilde conjugation rules and Eq.(§). 

Now we are going to explore these generators (A), as elements of a Lie algebra. Consider 
i = {cii,i = 1,2,3, ...} a Lie algebra over the (real) field R, characterized by the algebraic 
relations Oj Aa,- = C^ak, where G R are the structure constants and A is the Lie product 
(we are using the convention of sum over repeated indices). Taking the TFD Hilbert space, 
say H w , as the space of representations for £, we write 

[A i ,A j ]=iC*A h . (is) 



\A l ,A 3 \=iD%A k , (16) 



[A i ,A j ]=iE*A k ,. (17) 

That is to say, physically Eq.(^) represents the Lie symmetries characterized by the 
structure constants (by definition, the hat operators A are the symmetry generators, a 
class of dynamical observables) . Since the non-hat operators (A) can be identified with the 
observables, the other possible class of dynamical variables, then Eq. flTB]) shows the way that 
generators of symmetries determine dynamical changes in the observables. The Abelian (or 
non-Abelian) nature of the observables in regard to the measurement process is establised 
by Eq. ([17D (this interpretation has been used to write Eq . (|P7|) as [Aj,A,-] = and so to 
derive a classical TFD structure [|,0). 

In this procedure, we have split the twofold structure of Ti w , by introducing a representa- 



tion for Lie symmetry, defined by the (TFD) algebra given by Eqs.(|T5| ) -(|T7| ) , which we denote 
by *£ (regarding the close connection of this structure with the w*-algebra, see Ref. ||). In 
the next section we use this representation theory to analyse a spinorial representation for 
the *-Poincare group, *p. 

Closing this section, let us observe that for a gauge symmetry, a set of generators of 
the gauge group will be added to the space-time symmetry group; and, of course, these two 
set of generators commute with each other. However, the implications of gauge groups as 
additional quantum numbers, such as color, etc, will apear in defining the Hilbert space. 

III. POINCARE GROUP AND DENSITY-MATRIX FIELD EQUATIONS 

Here we build a *-p Lie algebra assuming that in Eqs.(|15|)-(|T7|) we have Df, = E\, = C^, 
where Ch are the structure constants of the Poincare Lie algebra, such that we can write 

[M^, P a ] = i(guaP„ - 9v»Pv), (18) 
[P„,Pu] = 0, (19) 
[M^, M ap ] = -i(g w M va - g up M^ + 9fi<T M pv - g vu M pix ), (20) 
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[M^, P a ] = [M fW , P a ] = i{g va P p - g aVL P v \ (21) 

[P^,Pu}=0, (22) 

[M^, M ap ] = -%{g pp M va - g vp M^ + g^Mpu - gu*M p(t ), (23) 

[M^, P a ] = i{g va P» - g ffli P v ), (24) 

[P„Pu] = 0, (25) 

[M pv , M ap ] = -i{g pp M uu - g up M pa + g^M pv - g uu M pp ), (26) 

where M pu stands for the generators of rotations and P p the generators of translations; 
whilst M pv and P p are, respectively, the corresponding observables for rotations and trans- 
lations. The metric tensor is such that diag(g pu ) = (1,-1,-1,-1), and g pu = for 
V v\ (i, v = 0, 1,2,3. 
The invariants of *p are 

W = Wpw", (27) 

P 2 = P^P", (28) 

W = w p w^ - w p w p = 2w p w p - w^vF, (29) 

p2 = (p 2 f = - p^pv = 2P p P tI - P P P P ; (30) 



where w p = \e pvpu M va P p are the Pauli-Lubanski matrices; e pvpa is the Levi-Civita symbol; 
and 

w p = X -e pvpa M va P p + \e pvpa M^P p ~ \e pvpa M VCJ P p ■ 
Notice that the tensor 

w p = ^e pupa M™P p 

can be used to define the scalar W = w p w p , which is not an invariant of *p but rather of 
the subalgebra (of *p ) given by Eqs.(j24|)-([26|). 

Denoting the vectors of the representation space of Ti w by \^f) such that 
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m = 10) ® (01 = 10) ® 10); (31) 

where |0) is an element of a Hilbert space H, such that Tt w = H®H*; we also assume that 
(0|0) = 1. The action of the operators O and O in Ti w is then 

0|*)->[O®l]|tf) 

= [O®1]|0)®(0| = O|0)® (0|, (32) 

and 

0|#)->[1® 

= [10 0] |0) (8) (0| = |0)® (010+. (33) 

As a consequence, we show that the invariant (P 2 ) has a null constant value in this 
representation, that is 

(P 2 ) = (2P^P M - P^l*) = 0. (34) 

This is but the density matrix equation for the Klein-Gordan field derived in /. 

In order to construct a spinor density matrix equation, we introduce an equation as 

KP M f = 0, (35) 

such that, differently from Eqs.(^), we write 

KP,f KP,f = (P 2 ). (36) 

Using Eq. (§), a generic solution is immediately found to be 

a" = <r<f, (37) 

where 7 M are the Dirac matrix, and a is some non null Lorentz invariant. A trivial choice is 
a = l (but another possibility for a is explored in the next section). Therefore we find 



( 7 "p„) |*) = ® i - 1 ® (yp^ 

= (7 M P M )|0) ® (0| - |0) ® (0|(7^) f = 0. (38) 
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where now > is a 16-component spinor, and |0)((0|) is the 4- (dual) Dirac spinor. 

Multiplying the rhs of Eq. fl38|) by \<p), it results in (•y^p^ — m) \<f)) = 0, the Dirac equation. 
Now, multiplying the lhs of Eq.(^) by (<f)\ it results in (0|(^7 M — m) — 0, the conjugated 
Dirac equation. In this sense, in fact, Eq.(^) is a Liouville-von Neumann-like equation for 
the Dirac Field. 

At this point, we are able to explore some symmetry properties of Eq. (|38D , as the 
invariance under similarity transformations. That is, consider 

( 7 ^)' = U (Y?,) U~\ (39) 

and 



\V)' = U\V). (40) 



Then Eq.(|38|) reads 



( 7 ^)' |*)' = 0. (41) 



On the other hand, if [U, (7 M -P M ) ] = then \^f)' given in Eq. (PH) is a solution of Eq. (P5|) . 
In this case, one suggestive example is provided by U = U{y^P^) written in the form 

U=U(yP)=exp[-iyP], (42) 

where y is the transformation parameter. Accordingly, Eq. (f40"D reads 

\V(y))' = exp{-iyP)\<f>) ® (0| exp(iyP t ). 

The gauge covariance can be considered if we write 

P -> -id ll = P ll + gA ll , (43) 

where is the usual covariant derivative. Then 

\V(y))' = exp(-yd)\<f>) <g> (0| exp^). (44) 



9 



where the operation ® is defined by Eg. (|31~D . This function \^(y)) is a solution of the density 
matrix equation, which is derived from Eg. (|38|) , but rather considering the gauge field, that 
is 

[(Yd,) <g> 1 - 1 <g> (Yd,)] |0) ® (0| = 0. (45) 

If we use the definition given by Eq . (j3l|) , then we can write (x,x'\^) = ^>(x,x'), such that 
Eq.(f44"l) can be written as 

ty(x,y) = exp [-yd(x)] (f)(x) <E> 4>(x) exp [yd(x)^]. (46) 

which is the generalized Heinz density operator. It is to be emphasized that \<f>) (or <p(x)) 
are elements of the Hilbert space. So the various quantum numbers either from the space- 
time group or from the gauge group are assumed to be attached to each element of the 
representation space. 

As another example of the use of such symmetry approach, we derive the TFD spinor- 
electro dynamics equations. The fields are assumed to be the basic dynamical variables, so 
that the Lagrangian may be written as 

L = 4>D(p + (p(D — D — m)$ 

+4>(D — D — m)<f> - 4>(D — D — m)4> 

2 v 4 l> ' 

which is derived from the usual spinor-electrodynamics Lagrangian and the hat mapping, 
and where 

D = %Y(d, + ieA,). 

Therefore, it follows that 

Dcf) + (D + m)4> -D4> = 0, 
(D + m)<j) = 0, 
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which are the quasi-particle TFD equations for the spinor field. The other equations can be 
derived accordingly. 

In the next section we are going to present another application of the method developed 
here by deriving a phase space Wigner structure attached to the Dirac field. We use the 
more general form of Eq.fl55|) including a ^ 1, in such a way that we can better compare 
our formalism with the one proposed by Bohm and Hiley |25| and Holland [2(| to derive 
phase-space spinor fields under a geometric perspective. 

IV. THE DIRAC FIELD AND A COLLISIONLESS BOLTZMANN EQUATION 

Notice that, using Eq. (p7|) , Eq. (|35|) can be written as 



(* l ^-v r Y^)*tsM = 0, (47) 

such the ^ts{x, x') = <f)(x)t4> s {x'), where the subscripts t, s = 1, 2, .., 4 refer to the two spinor 
indices. 

Here 7 ^ = 7^ ® 1, Y^ 1 = 1 ® 7 M are 16 x 16 matrices and 7 M are the Dirac matrices; 
a 1 = a ® 1 and o r = 1 C§> a; with a, the arbitrary Lorentz invariant to be specified; besides 
the 7 r,z -matrices fulfil two Clifford algebras: 

l lu } = {7 r/ \ l ru } = 2<T; {7" 1 , l rv } = (48) 
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Multiplying Eq.(f|7|) by (a 1 ) , it follows that 

A(k)% s (x,x') = (t^— - ^—^(xrf = 0, (49) 

where k = (a l )~ l a r . Using the following transformation 

did did 

77777 = ^(777^+^ (50) 



dx^ a/2 dq» M7 ' dx't 1 y/2 ' dq 
the square of A(k) can be written as 

A2 (-) = A r(-)^7 + + at(«) P ^ + Ar»P„P„ (si) 
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where 



A III/ / \ 

Af («) 


2 W 


VII \ / 1 u 

-«7 M )(7 


- K7 riy 


AT(«) 








AT(«) 


2 W 






AT(«) 


2 W 


+ «7 r ")(7 b ' 


+ «7 ri/ 



From these expressions for (/c), (i = 1, ..,4), the following operators can be defined 
o"+(k) = + «ry r ") a"(«) = - «ry"*); (52) 

then Eq.fl5l"l) is rewritten as 

A 2 (k) = ^(«K(«)^7 + a"-(KK(«) A 

Notice that if k anticommutes with the matrices j 1 ^ and 7 rAt , and satisfies n 2 = 1 (which 
is compatible with Eq.fl3"8|)), it leads to 

{a", a"} = {a" + ,a 1,+ } = 

Observe that the requirement of Lorentz invariance for k can be achieved if we define a 
in Eq. fl3"7f ) as a = 7 5 such that k = 7' 5 7 r5 , where 7 5 = Z7 X 7 2 7 3 7 4 . Using these results, we 
follow Bohm and Hiley to write 

A 2 = a^ + a u ^- p u + a»a» + Pfi ' 



pi*. 
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Therefore, from Eq. ([!£]), we obtain 



p^V(q,p) = 0; (53) 



12 



which is the Bohm and Hiley results | 25| , as developed by Holland p6| , about the relativistic 



phase space. But here, we have derived Eg. (|53|) from first principles, so that we can go 
somewhat further with it. Observe, in particular, that Eq.(^) can be derived also for the 
Klein-Gordon field if we use the transformation given by Eq.(|50|) in Eq. (|34]) . It should be 
stressed here that within an approach based entirely on representations of symmetry, only a 
collisionless Boltzmann equation can be obtained. The role of dynamics can be included by 
postulating the interaction Lagrangian, that will lead to collision terms. The result, then, 
will be model-dependant. Consistently, Eq. (|53"D describes only a drift term, that should be 
the same for any kind of systems, including the classical ones. 

In particular, in this formalism we take advantage of the fact that the state is treated as an 
amplitude, to implement higher order terms in Eq.(^) looking for an invariant Lagrangian 
procedure, avoiding then the (usual) intricate way to do that. Indeed, observe that one 
simple way to derive a Vlasov term (the influence of a external field via F^) in Eq. (|53|) , 
considering the first order (classical) limit in the Heinz approach, is to write the Lorentz 
invariant Lagrangian as 

-^*(q,p)mF,d^(q,p) - h(q,p)mF»d^*(q,p). (54) 

Then the energy-momentum tensor is defined as 

= [ d 3 p-y*(q,p)p^(q,p), (55) 
J p 

which is still compatible with the usual definition of T^ v since ty*(q,p)ty(q,p) can be inter- 
preted as the classical distribution function in the relativistic phase space. 

Another aspect to be explored in the present formalism is the inclusion of gravitational 
effects in the kinetic theory. The treatment of transport theory in a curved space-time 
background is somewhat hindered by the fact that a definition of the Wigner function 
ty(q,p) depends on the use of the Fourier transform of a spacetime correlation function 
with a translated argument. Fourier transformation (nor translations) are globally available 
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in general curved spacetime. However, as shown by Calzetta, Habib and Hu |23j , the 



construction is still meaningful, if carried out in the tangent space of an event in spacetime 
and related to the neighbourhood of that event by an exponential map. Therefore, starting 
from Eq. ([45]) , the Liouville-Vlasov equation for the Wigner function of a spinor field coupled 
to a gauge field with field strength tensor is then 

(- Pll F^D v +p»D ll )V(q,p) = 0, 
where _D M is the covariant derivative including the curved space-time spin connection. 

V. FINAL CONCLUSIONS AND REMARKS 

Summarizing, we have used here a representation formalism for Lie algebras, in which 
the elements of the representation space, say TC W , are closely connected with the notion 
of density matrix through the thermofield dynamics approach. We have started with the 
Poicare Lie algebra in the H^-space. Then we have derived a gauge invariant density matrix 
equation for the spinorial field, and shown that one particular solution for such an equation 
is the gauge invariant density operator as proposed by Heinz. This fact has set forth a solid 
structure to support the Heinz approach, since it has been derived here by an analysis of 
group theory, i.e. the fundamental symmetry properties of the system. In the same context, 
analysing a Lagrangian formalism, we have derived the TFD equations; and so, we have 
presented an association of the Heinz density method and TFD. The phase space structure 
of the spinorial field has also been studied as another application of our method. In this 
case, once again, we have been able to advocate the group theoretical bases to derive the 
early results of Bohm and Hiley [^5| and Holland ||2(J about the connection of the Clifford 



algebra and phase-space. But here we have proceeded somewhat further, pointing out a 
possible Lagrangian formalism to take into account higher order terms in a kinetic equation. 
This last fact opens the possibility to deal with the kinetic theory as a genuine quantum or 
classical field formalism (observe that, as we are analysing representations of a Lie algebra, 
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the function ty(p,q) in Eq.(|53D can be considered as a quantum-like field defined in a Fock 
space, similar to what has been developed for Galilean classical systems H0]); and so, a 
pertubative scheme can be introduced by symmetry properties of the Lagrangian, as has 
been argued in writing Eq. (|5^) . Besides, we have discussed the effect of the curved space- 
time background for the spinorial density matrix equation. Some of these aspects are treated 
in a more detailed way in another paper f27|). 
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